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Abstract. We show that fluid stationary models are able to reproduce the observed, negative vertical gradient of 
the rotation velocity of the extra-planar gas in spiral galaxies. We have constructed models based on the simple 
condition that the pressure of the medium does not depend on density alone (baroclinic instead of barotropic 
solutions: isodensity and isothermal surfaces do not coincide). As an illustration, we have successfully applied 
our method to reproduce the observed velocity gradient of the lagging gaseous halo of NGC 891. The fluid 
stationary models discussed here can describe a hot homogeneous medium as well as a "gas" made of discrete, 
cold Hi clouds with an isotropic velocity dispersion distribution. Although the method presented here generates 
a density and velocity field consistent with observational constraints, the stability of these configurations remains 
an open question. 
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1. Introduction 

Observations at various wavelengths show that some 
spiral galaxies are surrounded by a gaseous halo. 
This extr a-planar gas is multiphase: it is detect ed in 
Hi (e . g.. ISwaters. Sancisi. fc van der Hulsd Il997lh Ha 



Rand 2000 t iRossa et al.ll2004l). and X-r ay observations 



Wang et allf200lUStrickland et al.ll2004h ." n particular 
high-sensi tivity Hi observations of edge-on galaxies like 
NGC 891 llSwaters et al.ll997tlFraterriali et all2004a|) and 
UGC 7321 l|Matthews fe Woodll2003j) reveal neutral gas 
emission up to large dista nces from the plane (e .g. see the 
Hi map of NGC 891 in iFraternali et~a"fl [2005) and the 
presence of a negative vertical gradient in the gas rota- 
tional velocity A similar decrease of the rotational veloc- 
ity with distance from the plane has also been qbsery ed in 
the diffuse ion ized gas hal o of NGC 891 l|Randlll997(l and 
of NGC 5775 fe.andll2000h . 

The two major issues regarding the extra-planar gas 
are those of its origin and of its dynamical state. These are 
strictly related. For example, the halo gas could be t he re- 
sult of cosmological accretion (e.g., seelBinnevl2005[l. or o f 
a galactic fountain llShapfro fc Fieldll97fiHBregmarJl980lh 



or of both. Different structures and kinematics would be 
expected for these cases. Thus, the study of the origin and 
dynamics of the extra-planar gas not only is important in 
itself, but may also provide a new insight on the forma- 
tion and the structure of spiral galaxies. Furthermore, it 
may bring new evidence on the vertical distribution of 
dark matter. In this paper we focus on the problem of the 
dynamical state of the extra-planar gas. 

Two "extreme" types of models have been considered 
for the extra-planar gas and in particular for its decreasing 
rotational velocity: the ballistic and the fluid homogeneous 
models. 

Ballistic models describe the gas as an inhomoge- 
neous collection of clouds, subject only to the grav- 
itational potential of the galaxy: for example, in the 
galactic fountain model ionized gas is ejected from the 
galactic disk due to stellar winds and supernova ex- 
plosions, and t hen cools and falls back ballistically 
feregmanl Il980t . These models are able to explain ver- 
tical motions of the cold (Hi) and warm (Ha) gas 
components observed in several spir a l galaxies (e.g. 



Fraternali. Poste rior fc gancisi | |2fjf)4b[ |Booms ma et alJ 
2005|) . However. ICollins. Benjamin, fc Randl l|2002() have 
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tried a ballistic model to describe the extra-planar ionised 
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gas of NGC 891 and found problems in reproducing the 
observed kinematics. These authors suggest that the dis- 
crepancy could be solved by considering "the presence of 
drag between disk and halo, such as through magnetic ten- 
sion or viscous interactions between clouds. Alternatively, 
an outwardly directed pressure gradient could explain the 
gas kinematics" . 

We note here that the concept of pressure in the 
physics of the interstellar medium is a complex one, as 
there are several contributing sources of pressure, e.g. 
thermal, kinetic fo r turbulent), magnetic, cosmic ray and 
radiation (e.g., see lBoulares fc Coxlll99(l and references 
therein). For simplicity, in the following we will restrict 
our discussion to thermal and kinetic pressure. In fluid ho- 
mogeneous models, the extra-planar gas is described as a 
stationary rotating fluid without any motion along the ra- 
dial and vertical directions, with the galaxy gravitational 
field balanced by the pressure gradient and the centrifugal 
force. Until now, this approach has not been fully explored 
in all its possibili ties, and only a few attempts have been 
made (e.g., see iBeniaminlEcJo^l . Here we extend a pre- 
limi nary analysis of fl uid homogeneous stationary models 
l|Barnabe et al.ll2005|) to explain the vertical gradient of 
rotational velocity observed in the extra-planar gas. In 
particular, we show that models in which the gas pressure 
does not depend on the density alone (baroclinic configura- 
tions) are able to reproduce the observed vertical gradient. 
In addition, we show that baroclinic so l utions could pro- 
vide the drag invoked bv ICollins et alJ l|2002() . This sug- 
gests that a correct description of the extra-planar gas dy- 
namics may be found in "hybrid" ballistic-fluid stationary 
models. Finally, we address the question of the physical 
interpretation of our solutions and suggest, as an alter- 
native to the hypothesis of a hot homogeneous medium, 
the possibility of a "gas" of cold Hi clouds described by 
the stationary (fluid) Jeans equations, with the sustaining 
pressure given by a globally isotropic velocity dispersion 
tensor. 

The paper is organized as follows: in Sect. [2]and|3| we 
briefly introduce the baroclinic solutions and derive sim- 
ple and general rules for their construction. In Sect.^we 
construct a fluid homogeneous and stationary model for 
the galaxy NGC 891, and in Sect. we discuss the results 
from an astrophysical point of view. In the Appendices a 
simple, fully analytical model of a gas distribution with 
low rotation at high z is presented, together with the nu- 
merical code adopted for the case of NGC 891. 



2. The fluid approach 

2.1. The fluid equations 

We consider a gaseous axisymmetric system in permanent 
rotation, under the influence of an axisymmetric gravita- 
tional potential $ to t(ii, z): because of the axial symme- 
try, all the physical variables depend only on the cylindri- 



cal coordinates R and z. The stationary hydrodynamical 
equations for the gas are then 
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where p, P and fl denote the gas density, pressure and 
angular velocity, respectively; the gas rotational velocity 
is given by v v — flR, while Vr = v z = 0. Note that <I>tot 
represents the total gravitational potential, including the 
gas contribution. Later on, we will assume the condition 
that the gas is not self-gravitating and therefore $tot = 
where the galaxy gravitational potential $ is the sum of 
the dark halo and the stellar disk potentials. 

In standard applications, as for example the setting up 
of initial conditions for hydrodynamical simulations, the 
above equations are solved adopting a barotropic pressure 
distribution and neglecting the gas self-gravity. Thus, one 
fixes the gravitational potential $ and a specific function 
P(p), and integrates the first of Eqs. with the bound- 
ary p(R,0) or imposing p(R, oo) = 0. The angular veloc- 
ity fl is obtained from the second of Eqs. QJ. This leads 
to cylindrical rotation, i.e. f2 = £l(R). In fact, accor ding 
to th e Poincare-Wavre theorem ijLebovit d Il967t iTassoull 
1980), cylindrical rotation is equivalent to the fact that 
the acceleration field at the r.h.s. of Eqs. (TJJ can be ob- 
tained from an effective potential <& G ff (see Eq. jHJ), or that 
the gas density and pressure are stratified on $ e ff, and so 
the fluid is barotropic. 

Cylindrical rotation is in disagreement with the ob- 
served vertical gradient of the extra-planar gas rotation 
velocity, and this would seem to argue against the appli- 
cability of fluid stationary models. However, in the next 
Section we will show that it is possible to construct baro- 
clinic equilibrium solutions with a negative velocity gra- 
dient along z. Note that baroclinic solutions have been 
studied in the past for problems ranging from geophysics 
to t he theory of sun s pots and to ga l actic dynamics 1 (e.g., 
see lRosselandlll926t ITassoull ll98(A IWaxmanl Il978l and 
references therein). 

2.2. Baroclinic solutions 

For simplicity we restrict ourselves to non self-gravitating 
gas distributions, even though several results hold also 
in the self-gravitating case. We fix the total gravitational 
potential Q(R, z) for the galaxy but, at variance with the 
standard treatment, we prescribe a density distribution 
p(R, z) for the gas vanishing at infinity. (We do not ad- 
dress here the much more difficult problem of a consistent 

1 In particular, isotropic axisymmetric galaxy models can 
be interpreted as baroclinic fluid configurations, and show 

streaming velocities often decreasing wit h z (e.g., see 

lLanzoni fc Ciottil liool and ICiotti fc Bertir] l"2005l for simple 
examples). 
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assignement of non-zero pressure as a boundary condi- 
tion.). The first of Eqs. Q is integrated for the pressure 
as 

P(R,z) = J P^jdz', (2) 

where we also assume P(R, oo) = 0. In general, the ob- 
tained pressure P (and the corresponding temperature 
T = fmi p P/kp) can not be expressed as a function of 
p only, and so the system is baroclinic. Accordingly, the 
rotational velocity field 



v 2 JR,z) 



RdP 
~p"dR 



R 



d<S> 
dR 



(3) 



depends both on R and z. The major problem posed by 
the construction of baroclinic solutions is the fact that, 
for an arbitrary choice of p and the existence of phys- 
ically acceptable solutions (i.e. configurations for which 
1)%, > everywhere) is not guaranteed. In fact, due to the 
arbitrariness of the chosen density field, a negative radial 
pressure gradient in Eq. Q can be dominant for some val- 
ues of R and z. Thus, before addressing the specific case 
of the extra-planar gas in NGC 891, in the next Section 
we present a few general results of a mathematical nature 
that will be used as guidelines in the construction of physi- 
cally acceptable baroclinic solutions, while in Appendix lAl 
we present a fully analytical, baroclinic toy-model as an 
example of the procedure described above. 

3. Simple families of solutions 

The starting point of the following analysis is obtained by 
combining Eqs. (TJJ) and (0 and integrating by parts with 
the assumption that P and pd^/dR vanish for z = oo, 
thus obtaining the exact relation 



P% 
R 



dp 9$ 
dRd? 



dp <9$ 
d~z~'d~R 



dz' = C[p^}. 



(4) 



This "commutator-like" rel ation is not new (e.g., see 
iRosselan dlll92rl IWaxmaifl 978 and, in the context of stel- 
lar dynamics, Hunter 1977); here we note that the positiv- 
ity of the integrand in Eq. (0J is a sufficient (but not nec- 



3.1. Gas density distributions with a factor stratified 
on the effective potential 

We now elaborate in more detail the general results of 
Eqs. Q)-©. Let us consider the factorized gas density 
distribution 

p(R,z) =h(R,z)p e ($ eS ), (7) 
where h is a non negative function and 

$ off = $- / n 2 (R')R'dR' 

is the effective potential associated with the total potential 
$ and with a prescribed cylindrical rotation law £l(R); Rq 
is an arbitrary but fixed radius. We assume that p e in 
Eq. J7J) is a solution of the equation 



(8) 



VP = -pV$ ofr 



0) 



with assigned P = P(p). For = one obtains a hydro- 
static (and therefore barotropic) solution in the potential 
$ (that we indicate with p^), while for fl(R) ^ one has 
a cylindrical rotation (barotropic) solution. 

In practice, by adopting factorization (7J| one modifies 
a gas distribution with a cylindrical velocity field: this 
approach is of obvious interest because families of p e (§ c ff) 
can be easily constructed (see Appendix [BJ) . Substituting 
Eq. JJJ in Eq. |T3J gives 



pv% 
R 



m\R) 



dh <9$ dh d<f> 
dRdz 1 ~ dz'dR 
dp, 



p e dz' 



h^fjdz'. 

dz' 



(10) 



Equations (JSJ and IjlUI) prove the existence of phys- 
ically acceptable baroclinic solutions. For example, p = 
Ph($) +Pc($cff) leads to 



p c R 2 n 2 (R) 

ph + Pe 



(11) 



Another case of physically acceptable solutions is p = 
A(R)ph($>), with A{R) increasing and approaching a con- 



essary) condition to obtain vl > everywhere. Therefore, stant value for i? > 1. In this case v A v decreases reaching 



physically acceptable solutions are obtained if one assume 
a potential $ for which d$/dR > and d$/dz > 
(the usual situation) and a density distribution so that 
dp/dz < and dp/dR > 0. 

In addition, the bilinearity of C[p, <&] can be used to 
construct more complicate solutions starting from simple, 
physically acceptable "building block" configurations. In 
fact, from Eq. (jJJ it follows that the rotational velocity 
associated with p — p± + p2 is 



r Pl+ P2 

where pxv^^/R = C[pi, $] and P2V^ 2 /R 
if $ = $i + $ 2 , then 



u <p,2- 



(5) 

C[p2,$]. Also, 
(G) 



systemic (zero) velocity at infinity 2 . 

We now apply the method we have just described to 
a couple of more specific astrophysically relevant cases: a 
homeoidally stratified potential and a razor-thin uniform 
disk. For simplicity, we assume in Eq. J7J a hydrostatic 
density factor Ph(&)- 

3.1.1. Homeoidal potential 

Let <!>(£) be an homeoidally stratified potential with i 2 = 
R 2 + z 2 lq\ and < q<j> < 1; two well-known ex amples 
arc the|B innevl l)l98l[) logarithmic potential, and lEvans! 
l|l994l) spheroidal potentials. We write Eq. Q as 



p(R,z) = A(R)B(m)p h (<S>), 
2 E.g., the function A(R) = R/(l + R). 



(12) 
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where m 2 = R 2 + z 2 /q 2 and < q g < 1; A(R) and B(m) 
are positive functions. From Eq. (|1L)|) we have 



R 



k-^) A{R)R L 

A'(R) 



x B'(m) M zf 



p h (<S>)&(l)-dz' 
m I 



1% 



and so v 2 > if 



B(m) /9h ($)$'(£)-dz', 



$'(0 > 0, A'(R) > 0, B'(m) < 0, 



< <7$- 



(13) 



(14) 



Note that /?h does not enter in the sufficient condition 
(|14|) . The condition on Battenings (q g < q<t> ) requires that 
the gas density distribution must be stratified on homeoids 
which are flatter than the isopotential surfaces. This is, 
therefore, always satisfied for a flat gas distribution in 
the spherically symmetric monopole-dominated far field 
of any finite mass system (e.g. a stellar disk). 

3.1.2. Razor-thin uniform disk 

We now discuss the case of a razor-thin uniform disk. 
What are the conditions for having physically acceptable 
(v 2 , > 0) solutions near the disk? We explore this issue by 
assuming p = A(R)p h ($) in Eq. Q (or B = 1 in Eq. [TJ] ) 
and 



$ = 2ttG£ z. 

From Eq. O (or Eq. Q3]) 



A'{R) 
A(R) 



(15) 



(16) 



i.e. A'(R) > is the necessary and sufficient condition to 
have v 2 > in such a case. 

The physical reason for this condition is very simple. 
A gas distribution p(R, z) not stratified on $ as given by 
Eq. i(l"5|) must be rotating and according to Eq. its 
pressure must be radially increasing. From Eqs. l|15f) and 
J2J it follows that in the present case the pressure is pro- 
portional to the gas column density. This means that, in 
a vertical gravitational field, v 2 > whenever the column 
density is radially increasing. Note that this trend is con- 
sistent with the radial H I density distribution observed in 
several spiral galaxies (see lCavatte et alll994l) . Obviously, 
the case presented by more realistic disks (as the exponen- 
tial disk in Sect.EJ) will require the explicit construction of 
the whole equilibrium solution to check for the positivity 
of v 2 . 



4. Application to NGC 891 

We now apply the general results of Sect.|21to the modeling 
of the extra-planar gas of NGC 891 in an attempt to repro- 
duce its major features and in particular the decline of the 
rotational velocity with increasing z. The rotational ve- 
locity field resulting from this modeling will be compared 
with t he H I rotation curves derived by iFraternali et alJ 
1 20051) . 



4.1. The galaxy model 

We consider a very idealized mass model for NGC 891, 
consisting of three components: an exponential stellar 
disk, a dark matter halo with an asymptotically flat ro- 
tation curve, and a centrally peaked density distribution. 
This model provides the gravitational potential support- 
ing a (non self-gravitating) baroclinic gas distribution. Its 
parameters are fixed to reproduce the observed rotation 
curve in the galactic plane and to obey other observational 
constraints specified in the following. 

The surface density of the stellar disk is 



(17) 



[ ' 2nR 2 6 

where Md and Rd are the disk mass and scale-length, re- 
spectively. Its gravitational potential is given by 

Jo(kR) e~ k W 



GM, 



R d Jo (1 + fc 2 ) 3 / 2 



dk, 



(18) 



where Jo is the zeroth-order Bessel function of the first 
kind, and R = R/Rd, z = z/Rd (e.g. iBinnev fc Tremainel 
Il987ft . The observed rotation curve is reproduced by 
adding a two-component mass distribution whose poten- 
tial is 



2 

GM 



In R 



R z 



In 



(19) 



The first component represents a IBinnev! l|l98ll) loga- 
rithmic dark matter halo with asymptotic velocity vq and 
i?h = Rh/Rd', the corresponding potential belongs to the 
family of homeoidal potentials considered in Sect. 13.1.11 
The second component, needed to reproduce the steep ris- 
ing of the rotat ion curve in the innermost galaxy regions 
(see Fig.^l, is a lJaffd fl983) spherically symmetric density 
distribution of total mass Mq and scale-length ro- 

The adopted values for the disk scale-length and cen- 
tral surface brightness (Rd = 4.40 kpc. unn =21.4 mag 
arcsec -2 ) were taken from IShaw fc Gilmorel l(l989h . We 
find that the observed rotation curve of NGC 891 is well re- 
produced for vq = 138 km s _1 , Rh — Rd, 9* = 0.71, Mo = 
2.40 x 10 10 Af Q , r = 1.15 kpc, and M d = 7.7 x 1O 1O M . 
The resulting disk B-band mass-to-light ratio is ~ 3.5 (in 
solar units). 

4.2. The gas distribution and its rotational velocity 
field 

Having fixed the galaxy potential well, we now turn to 
the choice of the gas density distribution. Following the 
arguments presented in Sect.El we adopt the trial function 



p(R,z) = 



Po Ri (R + R) a 



,—*/h, 



R% (R 2 



(20) 



where m 2 = R 2 + z 2 lq 2 %1 R m = R m /R d , Rq = Ro/Rd, Po 
is the central gas, and (3 > a + 2 so that the total gas mass 
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Fig. 1. The adopted mass decomposition for NGC 891. 
Solid dots represent the observed Hi rotation curve 
( Fraternali et alJ 2005) , and the solid line the model rota- 
tion curve. The dashed, dotted, and dot-dashed lines show 
the separate contributions of the disk, the dark matter 
halo, and of the spherical Jaffe mass component, respec- 
tively. 

converges. Since the logarithmic halo belongs to the fam- 
ily of homeoidal potentials, p is a factorized distribution 
similar to that of Eq. H12[) with an homeoidally stratified 
component (flattening q g < a cylindrical radially in- 
creasing factor (a > 0) and an isothermal distribution 
stratified on the gravity field of a razor-thin uniform disk. 

We now solve Eqs. Q for the present model. In par- 
ticular, Eq. (J2J) is integrated numerically with a finite- 
difference scheme ( Appendix IHll . obtaining the gas pres- 
sure, velocity and temperature values on the numerical 
grid. 

The first result is that a large number of exploratory 
models (but not all of them), in which the value of the pa- 
rameters in Eq. (120ft is arbitrarily chosen, produce physi- 
cally acceptable rotational velocities. This confirms the va- 
lidity of the preparatory analysis in Sect. 01 For example, 
models without the radially increasing factor (obtained by 
setting a = in Eq. [20] ), invariably have < in the 
central regions, while models less flat than the halo po- 
tential (q g > g$) turn out to be non-physical in the far 
field. We conclude that the "core" properties of the phys- 
ically acceptable models are the gas radial depression in 
the central regions and a gas distribution flatter than the 
halo potential at large distances. The remaining parame- 
ters play only a minor role. 

We now present and discuss a specific, physically ac- 
ceptable model, which we refer to as the "Reference 
Model" (RM). The free parameters of the RM have been 
fixed to reproduce the o bserved Hi rotation curve of 
NGC 891 at z = 2.6 kpc JFraternali et alJl2f)05l) and the 
resulting rotation curve is shown in Fig. 01 (upper right 
panel, dotted line). We have obtained an acceptable curve 
by adopting a = 1, = 3.5, Rq — 1, R m — 2.25, h e = 1.5 
and q g — 0.1. Note that, due to the assumption of a neg- 



ligible self-gravity for the gas, from Eqs. J5J) and J2J it 
follows that v v and the temperature are independent of 
the specific choice of po- Other arguments will be used to 
fix the value of the latter (see Sect. I4.3)l . 

In Fig.OJwe show the radial profiles of the RM gas den- 
sity distribution at 5 different heights above the galactic 
plane (Fig. 0b) and the meridional sections of its isorota- 
tional (Fig. 0b) and isothermal (Fig. 0b) surfaces. ^From 
Fig. 0b it can be clearly seen that the rotation velocity 
decreases with z at a fixed R, while Fig. 0b shows that 
the gas is warm, with temperatures ranging from ~ 10 4 
to ~ 10 6 K. The hotter gas is near the R = axis, in 
correspondence with the density decrease near the galaxy 
center, while the lowest gas temperatures are attained near 
the galactic plane. 

In Fig. 01 we show the observed Hi rotation curves at 
four different heights above the galactic plane (z — 1.3, 
2.6, 3.9, and 5.2 kpc, dots), and the rotation curves from 
the RM (dotted lines), corresponding to horizontal sec- 
tions of Fig. Ob- The RM gas parameters have been tuned 
to obtain a good fit of the observed rotation curve at 
z = 2.6 kpc. This choice of parameters produces good fits 
for the curves at higher z too. It is remarkable that the 
gas density distribution in Eq. I|2U|) . mainly built on the- 
oretical arguments, should lead to a predicted rotational 
velocity decrease for the extra-planar gas which is so close 
to the observations. As to the discrepancy between data 
and model predictions at z = 1.3 kpc, it could be real (see 
Sect. l4.3|) . but it is possible that the observed data are af- 
fected by the limited angular resolution of the observations 
(HPBW = 28" w 1.3 kpc). Details on the derivation of the 
rotation curves from the observations will be presented in 
a forthcoming paper (Fraternali, in preparation). 

How do the RM properties depend on the specific 
choice of the gas density distribution? Among the param- 
eters, the largest effect is produced by the oblateness q g 
and the vertical scale-height h g . Increasing q g has a ma- 
jor effect: isorotation curves would be heavily modified for 
z > i?d/2, becoming almost vertical (cylindrical rotation) 
when q g ~ q$, which is the limiting flattening accord- 
ing to condition (14). We note that a similar behaviour 
is also shown by the analytical toy-model presented in 
Appendix E] even though this is based on a substantially 
simpler gas and galaxy model. A decrease of h g has a 
strong effect for z > Ra, producing a cylindrical rotation 
pattern extending out to R ~ 5i?d- 

Variations of a, /3 and Rq produce minor effects, at 
least in the regions covered by observations (R < 5Rd 
and z < 1.3i?d)- In particular, v v increases for increasing 
a and decreasing Ro, but the same overall structure of the 
RM is kept, and the differences are hardly detectable for 
z > 1.3i?d- An increase of j3 has an effect very similar to 
increasing a, and again there is a corresponding behaviour 
in the analytical toy-model of Appendix El Finally, an in- 
crease of g$ (i.e., the adoption of a rounder dark matter 
halo) at fixed q g has only a very marginal effect on v v . 
We conclude that baroclinic solutions are sensitive to the 
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Fig. 2. Extra-planar H I rot ation curves for NGC 891 (dots with error bars) at various heights over the galactic plane 
(from iFraternali et alJl200fit> . The solid curve is the rotation curve in the galactic plane (z = 0) as obtained from the 
model (see also Fig. The dotted lines are the predicted baroclinic extra-planar rotation curves of the Reference 
Model. 



adopted gas distributions, but only marginally to the halo 
flattening. It seems, therefore, that according to the fluid 
stationary interpretation presented here the kinematics of 
the extra-planar gas can not be used as a sensitive di- 
agnostic tool to investigate the flattening of dark matter 
haloes. 

4.3. Astrophysical interpretation 

The gas kinematics from the Reference Model is remark- 
ably similar to that observed for the extra-planar Hi; yet, 
the model gas temperature is much higher than that of 
Hi (Fig. 13;). Does this mean that fluid solutions should 
be abandoned as unphysical? 

The gas density distribution in Eq. H2()|) may be inter- 
preted not only as that of a smooth, homogeneous fluid 
(for which pressure and temperature are the usual ther- 
modynamical quantities), but also as the fluid description 
of the H i clouds distribution. In this paper we follow the 
first interpretation; however, a brief discussion of the sec- 
ond interpretation, in which Eqs. Q are interpreted as 
the Jeans equations for a system with isotropic velocity 



dispersion, is given in Sect. with further comments 

in Sect. 

4.3.1. A homogeneous gas distribution 

Equilibrium configurations (baroclinic or not) of a homo- 
geneous extra-planar gas are expected to be hotter than 
the Hi because of the needed vertical pressure support 
against the galaxy gravitational field. Therefore, the ob- 
served cloudy extra-planar H I gas can not be described di- 
rectly with fluid homogeneous models. We note, however, 
that baroclinic solutions as presented above are compara- 
ble to (or better than) ballistic models at reproducing the 
observed vertical gradients in the H I and Ha rotational ve- 
locities. Furthermore, there is, in addition to cold Hi, also 
hot gas in the halo of spiral galaxies a s observed in some 
edge on systems, including NGC 891 l|Bregman k, Pildisl 
ll994HStrickland et alJl200^ . It is therefore important to 
explore the possibility that the observed Hi traces the 
kinematics of an underlying, homogeneous, hot gas. 

The origin and the coexistence of hot and cold gas 
phases and their kinematical coupling are central issues 
in the discussion of possible realistic models for the extra- 



Barnabe et al.: Hydrostatic models for the rotation of extra-planar gas 



7 





z 

= / 












b J 


- / / 
: / / 

- / 








______ \60 

















































R H. 





R/R. 



R/R« 



Fig. 3. Panel a: radial profile of the density distribution 
for the Reference Model (RM) at z = 0, 1.3, 2.6, 3.9 and 
5.2 kpc (same heights as in Fig.[2J). Panels b,c: RM isoro- 
tation curves v v — const (in km s _1 ) and isothermal con- 
tours (in K) in the meridional plane. 



Fig. 4. Panel a: Bolometric surface brightness isophotes 
(in erg s _1 cm -2 ) for the RM. Panel b: meridional sections 
of equal cooling time (in yrs) surfaces for the RM with 
pa = 2.4 x 10~ 26 g cm~ 3 . Panel c: contours of equal £ = 
£cooi/£dyn (Eq. |27j ) in the meridional plane. 



planar gas. An important question in this respect is how 
long it would take for the hot gas to cool if not resup- 
plied with fresh energy, or to drag H I clouds (accreted or 
produced by galactic fountains) in regular motion around 
the galaxy axis. Clearly, hydrodynamical numerical simu- 
lations are needed to fully address it. Here, we limit our- 
selves to a first order analysis of the main astrophysical 
aspects and implications of the Reference Model presented 
above. 

Information on the thermal evolution of the RM is ob- 
tained from its local cooling time 



sity of electrons and protons, respectively. We also adopt 
the cooling function (in erg cm 3 s" 1 ) 



Ao(T) 



5.36 x 10" 27 T, 
2.18 x 10-18T- - 6826 



2.71 x iQ-^T 2 - 976 , 



(22) 



where the first equation holds for 10 4 < T < 1.3 x 
10 5 K, and the second for 1.3 x 1 5 < T < 10 8 K 
l|Mathews & Bregmanlll978l see also ICiotti et al.lll99l|) . 



For simplicity we assume n e 



E = ^A (T). 



(/i = 1/2) and so 



(23) 



tcoo\(R, z) = 



E 1 
— cx — . 

E Pa 



(21) 



where E = 3pkT/2pm p and E = n c n p A (T) are the gas 
internal energy and cooling rate per unit volume; k, p,, 
m p , n Q and n p are the Boltzmann constant, the gas mean 
molecular weight, the proton mass, and the number den- 



Note that, at variance with v v and T, E and i C ooi do 
depend on the value of the gas density normalization con- 
stant po in Eq. (|20|l . and so we are now forced to fix its 
value. A simple argument can be based on the require- 
ment of stationarity: the gas distribution must radiate per 
unit time a total amount of energy of the same order of 
magnitude as that provided by supernova explosions and 
stellar winds in the stellar disks of NGC 891, estimated 
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w 2.9 x 10 42 erg s" 1 (|Bregman fc Houcldll997l and ref- 
erences therein). The total (bolometric) luminosity of the 
gas distribution is 



47T 



ERdRdz cc Pq, 



(24) 



and from numerical integration we find that po w 2.4 x 
10 -26 g cm -3 , corresponding to a total gas mass M gas w 
1.6 x 1O 9 M (it should be noted that use of n c = 1.2n p 
to include helium ionization would lead to a density nor- 
malization and to a total mass w 10% smaller than the 
values obtained here). The hypothesis of a negligible gas 
self-gravity on galactic scale is thus satisfied, being the 
resulting gas mass only a few percent of the total mass 
of the galaxy. Curiously, the mass of the extra-planar gas 
distributed according to (HOJ (~ 5.5 x 10 8 M Q for z > 1.4 
kpc), is found to be of the same order of magnitude as 
the mass of the obs erved extra-planar Hi (w 6 x 10 8 Mq, 
ISwaters et al Jl997l) , while the luminosity of the RM extra- 
planar gas is very low, summing up to only « 5% of the 
estimated energy injection rate. From this point of view, 
the bulk of the supernova and stellar wind heating is ra- 
diated by the gas near the galactic disk, in good quantita- 
tive a greement with the conclusions of iRead fc Ponman I 
lll995h . This is confirmed by the surface brightness distri- 
bution E of the RM, obtained projecting E(R, z) along the 
line-of-sight. When the galaxy is seen edge-on, the surface 
brightness is given by: 



E(R,z) = 2 [ 

JR 



E R' dR' 



ir \jR a - R 2 ' 
In Fig. 0k we show the RM bolometric surface bright- 
ness distribution. We have also computed the RM lumi- 
nosity for the gas with T > 5.5 x 10 5 K (that we iden- 
tify with EUV-emitting gas), obtaining Leuv ~ 10~ 2 £ to t, 
with the bulk of the emission confined in a hot bubble 
along the galactic axis and matching the temperature dis- 
tribution in Fig. 

The conclusion that a major fraction of the galaxy 
heating is radiated near the disk is also confirmed by the 
fact that cooling times increase from the galactic plane to- 
ward regions of higher z and lower gas density, with values 
in the range 10 6 < t cool < 10 11 yrs (Fig.|U>). Obviously, 
the gas in the regions where t coo i is short will not re- 
main homogeneous and its natural fate will be to cool and 
form clouds and filaments. The energy input from the disk 
would be then used mainly to maintain a time-average, 
stationary multiphase structure. 

Additional insight on this multiphase medium is pro- 
vided by the dynamical time 

t dyB (R,z = 26 
v v (R, z) 

i.e., the orbital period of a gas element in circular orbit 
at (R, z). It is natural to introduce here the dimensionless 
number 

^cool 



1 

CX . 

£dyn Po 



(27) 



From Fig. it appears that £ > 1 for z > 1 -j- 2 R& 
and near the galaxy axis. It is expected that where £ < 1 
the gas will cool locally. Note that the regions where £ < 1 
match very closely those in which i coo i < 1 Gyr. 

In order to address the issue of the interaction of 
Hi clouds with the hot and homogeneous baroclinic gas 
distribution of the RM, we have computed the cloud 
drag time idrag, i-e. the time required for an Hi cloud 
to corotate with the homogeneous gas halo described by 
Eq. 1121)1) . This quantity is relevant whatever the origin 
of the Hi clouds: thermal instabilities of the hot gas, 
galactic fountains or cosmolo gical accretion. Following 
iKlein. McKee. fc Colellal l|l994j) . we define 



idrag(-R, z) 



8 r c 1 
— x oc — , 

oOd Vrcl PO 



(28) 



where Cd — 1 is a numerical coefficient, r c is the radius 
of a typical H I cloud, v rc i is the modulus of the relative 
velocity between the cloud and the homogeneous extra- 
planar component, and x = Pel P is the ratio between the 
cloud and the medium densities. 

We estimate tdrag assuming pressure equilibrium be- 
tween cold and hot gas components, from which x = 
T/T c rj 3000. A fiducial value of the relative velocity is 
obtained by assuming v rc \ = v v ~ 2 x 10 7 cm s _1 , while 
the cloud radius is estimated as 47rr 3 p c /3 = 10 5 MqM^, 
where M 5 is the cloud mass in units of 1O 5 M . We thus 



(25) obtain 



^dra 



8 / 3 x 1O 5 M M 5 
3v v \ 4im p m p 



1/3 



2/3 



2.7 x 10 s 



M 5 



1/3 



yrs, 



(29) 



where p = n p m p . From this simple formula it is appar- 
ent that massive clouds at high z above the galactic plane 
cannot be dragged by the hot and tenuous RM distribu- 
tion (see Fig. [3^ for characteristic values of n p ), while 
smaller clouds (e.g. M5 ~ 10~ 3 ) near the disk (z ps 1 -f- 2 
kpc) have drag times significantly shorter than rotational 
times. It cannot be excluded, therefore, that extraplanar 
H 1 clouds (or, even better, H 1 filaments, which have larger 
cross-sections) could be the tracers of the underlying fluid. 
Of course, if Hi clouds originated from cosmological ac- 
cretion, then the proper characteristic time to be com- 
pared with tdrag is the free-fall time. However, tdyn is of 
the same order of magnitude and the above arguments 
should still be valid. In any case, only hydrodynamical 
simulations will possibly clarify the importance of drag 
and related issues, such as the angular momentum redis- 
tribution between the homogeneous gas and the H 1 clouds. 
In conclusion, our analysis seems to leave open the pos- 
sibility advocated by Collins et al. of the presence of a 
non-gravitational forces acting on the extra-planar H 1 gas. 



Barnabe et al.: Hydrostatic models for the rotation of extra-planar gas 



9 



4.3.2. A "gas" of H I clouds 

As pointed out at the beginning of Sect. from a for- 
mal point of view, Eqs. JU are identical to the fluid 
stationary Jeans equations for an axisymmetric system 
with a globally isotropic velocity dispersion tensor. Thus, 
also for such a system we would obtain rotational ve- 
locities decreasing with increasing z. The pressure field, 
however, must be interpreted as P — per 2 , where a is 
the clouds (1-dimensional) velocity dispersion, so that 

cr = ^/kT/fim p w 41 \/f^ km s" 1 , where f 5 = T/10 5 K. 
From Fig-Efc it appears that near the disk the clouds veloc- 
ity dispersion would be of the order of 30 km s~ x or even 
less, while higher values of a are expected near the galaxy 
center. Such values might not be unreasonable. However, a 
further comparison of the hot and homogeneous gas solu- 
tion with the present one is not straightforward: while the 
RM gas distribution in Eq. (121 )[) was adopted to describe a 
hot and homogeneous gas distribution (and thus not very 
well constrained by observations), a proper Jeans-based 
analysis would require the choice of a cloud distribution 
similar to that observed. In principle this should be sim- 
ple to do, but we do not explore this any further here. 
At any rate, even if the temperature and related problems 
addressed in Sect. 14.3.11 were solved, in the Jeans-based 
interpretation the question of how the clouds velocity dis- 
persion tensor is rendered isotropic would become the new 
central question. 

5. Discussion and conclusions 

We have presented a family of fluid stationary models for 
the extra-planar gas in spiral galaxies. As an application, 
we have built a model for the extra-planar gas of the spiral 
galaxy NGC 891. The main results of our analysis are: 

(i) Physically acceptable baroclinic solutions exist for gas 
density distributions (isodensity surfaces) more flat- 
tened than the isopotential surfaces. For regions near 
the disk, the condition is that the gas density distri- 
bution is centrally depressed. 

(ii) Application of our method to the case of the edge-on 
galaxy NGC 891 has shown that with baroclinic solu- 
tions it is possible to reproduce the observed vertical 
decrease of the gas rotational velocity. 

(iii) In homogeneous fluid stationary solutions the gas tem- 
perature is in the range 10 4 < T < 10 6 K, well above 
the Hi temperature and the cooling times are short, 
of the order of the orbital times. 

(iv) In the hot, homogeneous gas configurations with nega- 
tive velocity gradients of the present models, the drag 
on small Hi clouds (~ 10 2 M Q ) near the disk is impor- 
tant and may account for the observed vertical velocity 
decrease. It is not certain, however, that this applies 
to the case of NGC 891 where the Hi structures may 
be more massive. 

(v) Instead of a hot homogeneous medium, the baroclinic 
solution can also provide a good fluid description of 



a "gas" of clouds (according to the stationary Jeans 
equations) with a globally isotropic velocity dispersion 
tensor. All the above considerations about the decrease 
of rotational velocity with increasing z would still hold, 
while the temperature field would have to be inter- 
preted as the velocity dispersion field. 

Ballistic models have been tried for the extra-planar 
gas. It seems, however, that in order to reproduce the ob- 
served rotational velocity gradients, also non-gravitational 
effect s, such as gas and magnetic pressure, may be neces- 
sary ijCollins et al.ll2002HFraternali k Binnev 112005^) . Our 
analysis seems to suggest that such non gravitational ef- 
fects could be due (at least for low mass clouds) to ram 
pressure. Note that £d rag depends on the relative velocity 
between the Hi clouds and the hot baroclinic gas, hence 
the effect of the drag would be to regularize the motion of 
the clouds to the velocity field of the embedding medium. 
For example, if the H I clouds are ejected from the disk by 
galactic fountains, they will be, on average, initially faster 
with respect to the baroclinic gas and will be slowed down 
to its velocity while rising. Incidentally, the drag could 
also be effective in regularizing the dynamics of accreted 
clouds: from this point of view, a regular Hi dynamics (in 
particular, the same sense of rotation of the galactic disk) 
would not be necessarily an indication of internal origin 
for the extra-planar gas. 

The problem of stability of the models presented here 
is a difficult one and remains open. Numerical hydrody- 
namical simulations would be very useful to address it. 
In the approach adopted in this paper, where the density 
distribution for the extraplanar gas is assigned, pressure, 
temperature and rotation cannot be arbitrarily prescribed 
but are determined by the galaxy gravitational potential. 
The stability of such configurations is not guaranteed. 
Additional effects worth studying and needing quantita- 
tive estimates are the drag of clouds (both of internal and 
external origin), their cooling and evaporation times, and 
the overall energy budget , to fo llow up on the pioneer- 
ing wo rk of lCox fc Smith I l|l974h and lMcKee fc Ostrikerl 
Equally important would be the numerical study 
of baroclinic solutions in the context of Jeans equations. In 
this case, isotropy of the velocity dispersion tensor could 
be provided by cloud-cloud collisions, a physical ingredi- 
ent usually neglect ed in numerical simulations (however, 
see IWaxmanlll978h . One could speculate that the prob- 
lems encountered by purely ballistic models is due to the 
absence of a collision term in the simulations, as already 
suggested bv lCollins etal] l|2002|) . 

We conclude remarking that fluid-homogeneous and 
Jeans-based baroclinic models do not exclude each other. 
Perhaps, they could be integrated and used together for a 
better understanding of the dynamics of the extra-planar 
gas. 
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Appendix A: A simple analytical toy-model 

We present here a simple and fully analytical toy- model of 
a baroclinic gas distribution characterized by the vertical 
decline of the rotational velocity. The dimensionless gas 
density distribution 



f2+A 



(0 < A < 3), 



(A.l) 



7 - ?7h(4 - 7) (7 - 1) 



VhR 2 



(5-7)(3-7)(7-2)fT- 2 (5-7)fT 



(A.2) 



(Ciotti & Bertin, Eq. [27]). The pressure (normalized to 
AnGRgPohPo), and the square of the gas rotational ve- 
locity field (normalized to 4irGR 2 poh) are obtained from 
Eqs. For simplicity we report here the solutions 

up to the linear terms in the flattenings, even though 
the derivation of the full solutions (i.e., including also the 
terms in ?7h?7g) presents no difficulty. Accordingly, 



P 



(5- 7 )(3-7) 



7 + A 



(3~ 7 )7^ + (5-7)%A i? 2 
7 + A f 2 



2XR 2 



(7 + A)fT+ A V3-7 



7 



Note that, for rj & < ?7h(3 — 7)/ (5 — 7) the (linearized) 
solution is unphysical. The isorotation curves for a physi- 
cally acceptable case (obtained from the full solution) are 
shown in Fig lA.ll where it is apparent how, even in this ex- 
tremely simplified model, the rotational velocity decreases 
with increasing z. 



15 



10 



is obtained from a homeoidal expansion of the oblate 
power-law distribution po/m x , where p = p/po, r = 
VR 2 + z 2 : m 2 = R 2 + z 2 /(l-r] g ) 2 , R = R/R g , z = z/R s , 
and < rip < 1/A is the flattening of the distribution 
(ICiotti fc Bertir]l2005h . 

The galaxy (stars plus dark matter) density distri- 
bution is also described by the homeoidal expansion of 
PoYi/fn 1 ', (where < 7 < 3, poh is a normalization den- 
sity, m is defined like above but with rfo replacing r) g , 
and < ?7h < 1/7)- Without loss of generality also the 
galaxy scale-length is R g , and the gravitational potential 
for 7 7^ 2 (in units of AnGR 2 poh) is given by 



(5-7)(l-%A)-(4-7)( 7 -l)7? h 



(A.3) 



(A.4) 
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Fig. A.l. Meridional section of the isorotation surfaces w 2 
(labels along the curves) for 7 = 2.5, % — 0.2, A = 2.3 
and r] g = 0.3. 

Appendix B: How to construct p c (& c s) 

A first approach to the construction of the distribution 
Pc^cff ) considered in Sect. 13. H is the assumption that P = 
P(p), as for example P cx p 1 . In this case, from Eq. 10 



Pe($eff) = POX < 



exp 



1 



pm p 
kT 

7 - 1 P m p 
7 kT 



($ eff - $ ff,o) 

($eff,0 - $cff) 



(B.l) 



where po and ^cff.o are taken at the same arbitrary but 
fixed point (Ro, Zo), and the first expression holds for the 
isothermal case. In the second case T/Tq = (Pc/po) 7-1 
and, at variance with the isothermal case, a truncation of 
p e may appear. A different approach to the construction 
of p c is however possible, where the specific density field 
is prescribed. This can be done by using the well known 
property that solutions of Eq. J^J necessarily are stratified 
on $ cff , i.e. p c = p c ($eff) and P = P($ of f). Thus, if one 
fixes $ c ff and a prescribed function p c (<I> c ff), the pressure 
field is obtained by direct integration as 



P($ cff ) = P($eff,0) 



Pe{t) dt. 



(B.2) 



Of course, this approach can be used only for density 
stratifications such that P > everywhere. 

Appendix C: The numerical code 

We describe here the scheme on which the (double pre- 
cision) FORTRAN77 code adopted for the construction 
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of models in Sect. 0] is based. From symmetry arguments 
we restrict the computation to the region [P m i n , P max ] x 
[z min , z max ] in the half-space z > 0, on which we define a 
bilogarithmic grid of (jir + 1) x (n z + 1) elements. In the 
radial direction the grid is 

Ri = Rix 10 {t - 1)AR , (C.l) 

where Po = 0, Pi — R min , and 
log if! 

max 

log P min 



AR = 



(C.2) 



n R - 1 

so that i? max = Rn R the z coordinate is discretized in the 
same way. 

The gas density and the gravitational field components 
at {R u z 3 ) are p), R $) = (d$/dR)ij and 2 $j = (dQ/dz)^. 
For fixed Ri we integrate the first of Eqs. Q with linear 
interpolation and boundary condition P(Ri, z max ) = 0: 



^ n z — 1 

P) = g £ x MU*fc + Pfc+i***+i)> 

fe=J 



(C.3) 



where P| = P(Ri, Zj) and Az^ = z^+i — Zk- Accordingly, 
from Eq. ©: 



— fni+l 



where R P] = (Pj +i - P))/Ai^ and AP; = 



(C.4) 
Pj+i — Pi- 

The gas temperature T, the cooling rate E (Eq. the 
cooling time t coo \ (Eq. [21| ) and the ratio £ between cooling 
time and dynamical time (Eq. |27| ) are then obtained on 
the grid points. 

The gas total luminosity is obtained as 



riR — l Tiz — 1 



I<tot - 51 

where 

rRi+i rZj+i 

L) = 4tt / / RE(R, z) dPdz, 



(C.5) 



(C.6) 



and E(R, z) is evaluated on the numerical grid from 
Eq. I|23|l . The integral in Eq. (|C.6|1 is evaluated analiti- 
cally on each region [P,-, P,+i] x [zj, Zj+i] by considering 
the bilinear expansion 



E(R, z) = A + AiR + A 2 z + A 3 Rz, 
where 

• ■ - E ) E )+i - E ) 
Aq = E\ — 3 A „ J P a - -2±± 3 - Zj + 



(C.7) 



AP 4 



AR,A 



z. 



A _ E J ~ E 3 + 1 - E J + E J - ^J + l 9 



ARiAzi 



Ao 



A, 



f; ' - ^ +1 +E) -E) +l 



Azj ARiAzj 

+ 1 ~ E 3 + % ~ E 3+l 

ARlAzh 



Ri 



With this formula, 
L) = TtARiAzj 2A Q (R l+l 



Ri) 



4Ai 2 2 

-^-(-Ri+i + Ri+iRi + Ri ) 

A2{Ri+i + Ri)(zj + i 
2A 3 
3 



(Pf +1 + P J+ iP 4 + P?)(^+i + 



(C.8) 



Finally, the edge-on surface brightness £ is given by the 
discretization of Eq. (|25|) and so, in our scheme, 



71R — X 



k—i 



Rl +1 -RI- JRI R? 



.(C.9) 
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